Isotope effect in quasi-two-dimensional metal-organic antiferromagnets 
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Although the isotope effect in superconducting materials is well- documented, changes in the 
magnetic properties of antiferromagnets due to isotopic substitution are seldom discussed and re- 
main poorly understood. This is perhaps surprising given the possible link between the quasi-two- 
dimensional (Q2D) antiferromagnetic and superconducting phases of the layered cuprates. Here we 
report the experimental observation of shifts in the Neel temperature and critical magnetic fields 
(ATn/Tn ~ 4%; AB C /B C « 4%) in a Q2D organic molecular antiferromagnets on substitution of 
hydrogen for deuterium. These compounds are characterized by strong hydrogen bonds through 
which the dominant superexchange is mediated. We evaluate how the in-plane and inter-plane ex- 
change energies evolve as the hydrogens on different ligands are substituted, and suggest a possible 
mechanism for this effect in terms of the relative exchange efficiency of hydrogen and deuterium 
bonds. 



In 1950 the isotope effect was first observed in the su- 
perconducting critical temperature of mercury [TJ [2] . The 
effect was key to verifying the electron-phonon pairing 
mechanism responsible for superconductivity in conven- 
tional superconductors and played an important role in 
motivating the microscopic theory of Bardeen, Cooper 
and Schrieffer [3 j. The analogous pairing mechanism in 
the superconducting cuprates is still open to vigorous 
debate. It seems increasingly likely that some kind of 
magnetic interaction is responsible, arising from the un- 
derlying quasi-two-dimensional (Q2D) Heisenberg anti- 
ferromagnetic lattice [U [5j El [3 El |9] . Nevertheless, an 
electron-phonon coupling component cannot yet be ruled 
out [TUJ [TTJ [12] . One reason for this is that an adequate 
explanation of the anomalous, doping-dependent isotope 
effect observed in the cuprates is not forthcoming p~3j [14] . 
Given these considerations the importance of achieving 
an understanding of the effects of isotopic substitution 
on the properties of magnetic systems of reduced dimen- 
sionality is clear. 

Metal-organic coordination complexes (see for exam- 
ple p~5j [16] [17] ) are an ideal framework in which to make 
studies of Q2D magnetism, and hence test theories relat- 
ing to unconventional superconductivity. The extremely 
anisotropic Heisenberg exchange interactions observed in 
these materials are highly susceptible to small perturba- 
tions [18 . Furthermore, the versatility of coordination 
chemistry allows for a high level of control over their 
composition and crystal structure. In consequence, the 
dimensionality of their magnetic structures can be readily 
changed via adjustments of molecular architecture and 
the magnitude of the exchange energies altered so that 
they are accessible by laboratory magnetic fields at cryo- 
genic temperatures [18]. Isotopic substitution provides 
an additional fine-scale tuning of magnetic properties. 
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FIG. 1: Structure of CuF2(pyz)(H20)2 at room temperature 
viewed (a) perpendicular to the a-axis, showing the Cu-pyz- 
Cu chains; and (b) parallel to a-axis, showing the 2D magnetic 
planes and the F-H-0 hydrogen bonds that mediate the dom- 
inant exchange interactions. Cu = red; F = green; N = blue; 
C = grey; H = cyan; and O = magenta. Hydrogen bonds are 
shown as dotted lines. 



We have chosen to measure the coordination com- 
plex CuF2(pyz)(H20)2 (pyz = C4H4N2) because its 
quasi-two-dimensional antiferromagnetic properties are 
directly analogous to the parent compounds of the 
cuprate superconductors, which are similarly based on 
square lattices of Cu 2+ ions. We observe for the first 
time in such a system measurable isotope shifts in both 
the Neel temperature and the critical magnetic field. 

The structure of CuF 2 (pyz)(H 2 0)2 and its isotopes 
(see Fig. [T]) is similar to the chain-like compound 
Cu(N03)2(pyz) [19] . In that material electron para- 
magnetic resonance (EPR) suggests that the occupied 
d x 2_ y 2 orbitals are arranged such that one of the lobes lies 
along the Cu-pyz-Cu chain direction [20]. This, together 
with a lack of a strong exchange pathway between the 
chains leads to Cu(N03)2(pyz) exhibiting extremely one- 
dimensional (ID) magnetic characteristics [20]. The situ- 
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ation in CuF 2 (pyz)(H 2 0)2 is very different; the presence 
of the strongly electronegative F and O ions cause the 
d x 2_ y 2 orbitals to lie in the O-Cu-F planes diminishing 
the exchange along the chains [19] • In addition, strong F- 
H-0 hydrogen bonds develop between the chains that are 
able to efficiently mediate superexchange. The net result 
is that the magnetic properties of CuF 2 (pyz)(H 2 0)2 ap- 
proximate to a quasi- two-dimensional spin-^ Heisenberg 
antiferromagnet on a square lattice [T9] . 

In this paper we show that isotopic substitution of hy- 
drogen for deuterium in CuF 2 (pyz)(H 2 0)2 has a sizable 
effect on its magnetic properties. In particular, deuter- 
ation of the water molecules acts to suppress the in- 
plane exchange energy without significantly changing the 
inter-plane exchange energy, the crystalline-electric-field 
local to the magnetic ions, or the crystal structure. This 
implies that hydrogen bonds are more efficient superex- 
change pathways than deuterium bonds at low tempera- 
tures. The lower mass leads to a higher zero-point vibra- 
tional energy of the hydrogen ion in the F- • • H — O bond 
compared to that for the deuterium in the F- • • D — O 
bond [21j [22] . We propose that this is responsible for the 
observed changes in antiferromagnetic ordering temper- 
ature and critical magnetic field. 

Single crystal samples of CuF 2 (pyz)(H 2 0) 2 were pre- 
pared via the recipe described in Ref [19] . Four different 
hydrogen/deuterium isotopes were produced: one con- 
taining only hydrogen; one in which the hydrogen ions 
on the pyrazine ligands are replaced with deuterium; one 
in which the water is replaced with D 2 0; and another 
containing only deuterium. For brevity the isotopes will 
henceforth be referred to as hg, aVpyz, d 2 o and dg : re- 
spectively. The substitution was achieved by using com- 
mercial pyrazine-a^ (Aldrich, 98 atomic % D) and D 2 
(~ 95 atomic % D). DC magnetization measurements 
were performed using a Quantum Design SQUID mag- 
netometer with the magnetic field B applied parallel to 
the a-axis. Pulsed- field magnetization experiments were 
performed using a technique described elsewhere [23] [24] . 
Fields were provided by the 60 T short-pulse magnets at 
the Nicholas Kurti Magnetic Field Laboratory, Oxford 
and the 65 T short-pulse magnets at National High Mag- 
netic Field Labratory, Los Alamos. The susceptometer 
was placed within a 3 He cryostat providing temperature 
T down to 500 mK. 

Typical low-temperature, pulsed- field magnetization, 
M(B), data are shown in Fig. [2] for the four different 
isotopes with the magnetic field applied parallel to the 
a-axis. The form of the data is characteristic of quasi-2D 
magnets (cf. the measured and simulated magnetization 
curves reported in Ref. [H]): a concave shape at low 
fields followed by a transition to saturation at critical 
field B c . The concavity of the data in Fig. [2] implies that 
these compounds have a very high degree of exchange 
anisotropy. Taking the critical field to be the position 
of the peak in the curvature of M(B) we find an aver- 
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FIG. 2: Pulsed- field magnetization M data versus magnetic 
field B for the four different hydrogen/deuterium isotopes of 
CuF 2 (pyz)(H 2 0)2. The data are taken at T = 1.5 K with 
B applied parallel to the a-axis and are normalized by the 
saturation magnetization M sat . hg = CuF2(/i4-pyz)(H20)2, 
d 4 -pyz = CuF 2 (^4-pyz)(H 2 0)2, d 2 o = CuF 2 (/i4-pyz)(D 2 0)2, 
and d 8 = CuF 2 (d4-pyz)(D 2 0)2. 



age B c for each isotope using a number of datasets from 
several different samples at T < 1.5 K. The values of B c 
extracted in this way are shown in Fig. [3|a). It is clear 
from the curves and the B c data that deuteration of the 
pyrazine rings alone has very little effect on the magne- 
tization. In contrast, deuteration of the water molecules 
significantly reduces the critical field. This suggests that 
changes in the in-plane exchange are largely responsible 
for the differences between the isotopes. 

A simple quasi-2D Heisenberg Hamiltonian [18] sug- 
gests that at the critical field g/iB^c = 4 J + 2J± [25] . 
Because of the large exchange anisotropy (i.e. J ^> J±) 
we can obtain an estimate for the in-plane exchange pa- 
rameter from the M(B) data. The results are shown 
in Fig. |3|d) and make use of the ^-factor obtained from 
electron paramagnetic resonance experiments on the hg 
isotope (g a = 2.42 ± 0.02 [H]). These estimates assume 
that the differences in B c between the isotopes are due 
almost entirely to differences in the in-plane exchange pa- 
rameter and that isotopic substitution of hydrogen has no 
significant effect on either the ^-factor or the inter-plane 
exchange energy. To check the validity of these assump- 
tions we make use of the low-field, a-axis susceptibility 
x(T) data shown in Fig. [IJ 

It is seen from the figure that the susceptibilities are 
broadly similar for all the compounds. The antiferromag- 
netic ordering temperature Tn is the sharp peak observed 
in x(T) at low temperatures and its value for the four iso- 
topes are plotted in Fig. [3jb), where it is seen that they 
follow the same trend as the critical field data. The value 
of Tn agrees well with the results of other techniques [19] . 

The xCO data are fitted for 5 < T < 50 K using the 
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FIG. 3: Isotopic variation in: (a) critical magnetic field; (b) 
antiferromagnetic ordering temperature; (c) the g- factor par- 
allel to the a-axis deduced from fits to x(T) data; (d) in-plane 
exchange energy deduced from fits to x(T) data (solid circles), 
and from critical field data (open circles); and (e) inter-plane 
exchange energy J± deduced from the expression in Ref . [27] . 



technique outlined in Ref. [26], which uses the results 
of theoretical predictions of the paramagnetic suscepti- 
bility of a quasi-2D spin-^ Heisenberg antiferromagnet 
with a single-exchange Hamiltonian. The fits yield values 
of J and g a . Analysis of quantum Monte-Carlo calcula- 
tions predicts that the inter-plane exchange energy J± for 
quasi-2D Heisenberg magnets is related to J and Tn via 
the expression |Jj_| = |J|exp(2.43 - 2.30 x \J\/T N ) [27]. 
Thus, using the parameters obtained from the x(T) data, 
the magnitudes of J±_ can be deduced for all the isotopes. 
The values of g a , J and J± are plotted in Fig. |3jc), (d) 
and (e). 

As expected from the form of the M(B) curves the 
values of J± are indeed small compared to J, leading 
to anisotropics of the order of 1 x 10 -3 . The ^-factor 
for the hS compound (g a = 2.40 ± 0.02) is in excellent 
agreement with that obtained from EPR experiments 
(#epr = 2.42 ± 0.02) [IS]. Although the fitted values 
of g a vary from isotope to isotope this variation is small 
(Ag a /g a « 2% on substitution of all the hydrogens for 
deuterium) and on a par with the experimental error. In 
addition, the small values of J± mean that the errors on 
Tn, J and g a , once propagated through the above expres- 
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FIG. 4: (a) Low- field susceptibility of the four different hydro- 
gen/deuterium isotopes of CuF2(pyz)(H2 0)2 measured with 
B = 0.05 T applied parallel to the a-axis (points) compared 
to the fits to the expression from Ref [26 (red lines), (b) Ex- 
panded view close to the ordering temperature Tn. The data 
are offset for clarity. 



sion, encompass its variation across the isotopic family. 
The results suggest that substitution of hydrogen for deu- 
terium in CuF2(pyz)(H20)2 has very little effect on the 
^-factors or the inter-plane exchange energy. Therefore, 
the observed shifts in the Neel temperature and the crit- 
ical magnetic field reflect the changes in the in-plane ex- 
change energy. The variation in J between isotopes with 
hydrogenated and deuterated water molecules is signif- 
icant (A J I J « 4%), whereas little or no difference is 
observed on deuteration of the pyrazine molecules alone, 
confirming the picture of the dominant superexchange 
being mediated by the hydrogen/deuterium bonds be- 
tween the carbon and fluorine ions. 

A shift in critical temperature is also ob- 
served on deuteration of certain layered organic 
molecular superconductors. Here, substitution 

of hydrogens on the ethylene end-groups of the 
bis(ethylenedithio)tetrathiafulvalene, or ET, molecule 
leads to an increase in the superconducting critical 
temperatures [28j [29]. This is caused by a shortening 
of the C-D bond compared to the C-H bond of the 
hydrogenated sample, so that the spatial extent of the 
conducting layers is decreased and the interlayer transfer 
integral is slightly suppressed [29j [30]. Exactly how 
this promotes the onset of superconductivity is not yet 
entirely clear, but it is suspected that reducing the 
dimensionality of the Fermi surface in this way enhances 
its susceptibility to the antiferromagnetic fluctuations 
believed to play a vital role in the formation of the 
superconducting condensate [30] . 

CuF2(pyz)(H2 0)2 is non- metallic and so Fermi sur- 
face effects have no part in the observed isotope ef- 
fect. Our preliminary neutron scattering data show that 
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there is very little difference between the hydrogenated 
and deuterated structures of this material. In partic- 
ular, the length of the F- • • H — O hydrogen bonds and 
the F- • • D — O deuterium bonds are identical (dp...o = 
2.61±0.01 A at room temperature). However, it is known 
that at low temperatures the effective "size" of deuterium 
in a chemical bond is smaller than that of hydrogen: the 
enhanced mass of the deuterium ion leads to a smaller 
zero-point vibrational amplitude compared to the equiv- 
alent hydrogenated bond [2TJ [22] . This localization leads 
to a reduction in the overlap of the electronic wavefunc- 
tions, a suppression of the transfer (or hopping) integral 
through the deuterium bond, and hence a decrease in the 
efficiency of the superexchange along this pathway. We 
suggest that this accounts for the observed shifts in Tn 
and B c in the isotopes of CuF2(pyz)(H20)2- 

This mechanism is similar to that proposed 
to occur in the hydrogen/deuterium isotopes of 
the quasi- ID Haldane-gapped antiferromagnet 
Ni(C 5 Hi4N2)2N 3 (PF 6 ) [31]. In the latter case, deutera- 
tion is believed to slightly diminish the efficiency of the 
secondary exchange path leading to changes in the mag- 
netic field/temperature phase diagram. Furthermore, a 
reduction in the Neel temperature (ATn/Tn ~ 0.6%) 
observed in La2Cu04 caused by substitution of O 18 
for O 16 is attributed to suppression of the zero-point 
motion of the apical oxygen ions [32]. However, due to 
the magnitude of J, it is not possible to measure B c in 
that compound. In contrast to those two materials, the 
isotope effect we have measured in CuF2(pyz)(H2 0)2 
is large and seen in both Tn and B c (ATn/Tn ~ 4%; 
AB C /B C « 4% for hS and d8). Nevertheless, the 
similarity between the proposed mechanism suggests it 
may be a common effect in low-dimensional Heisenberg 
antiferromagnets. Exploitation of the effect could allow 
materials engineers to fine-tune exchange interactions in 
solid-state magnetic devices. 

Finally, we note that a recent study has described how 
the effective sizes of bonded hydrogen and deuterium in- 
crease with temperature, and at different rates [21J. The 
implication is that for some materials, whose magnetic 
properties depend crucially on superexchange through 
hydrogen bonds, it may be possible to observe an ex- 
change energy that is enhanced at higher temperatures 
or an isotope effect whose sign changes as a function of 
temperature. 

In summary, we have measured significant shifts in 
both the ordering temperature and critical magnetic field 
in a quasi-two-dimensional Heisenberg antiferromagnet 
due to isotopic substitution. We associate these shifts 
with differences in superexchange efficiency between hy- 
drogen and deuterium bonds. The results highlight the 
fact that relatively small changes in composition can give 
rise to large perturbations in the properties of highly 
anisotropic magnetic materials; an important consider- 
ation when attempting to understand the anomalous iso- 



tope effect in the cuprate superconductors. 
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